Pauli paramagnetism of an ideal Fermi gas 
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We show how to use trapped ultracold atoms to measure the magnetic susceptibility of a two- 
component Fermi gas. The method is illustrated for a non-interacting gas of 6 Li, using the tunability 
of interactions around a wide Feshbach resonances. The susceptibility versus effective magnetic 
field is directly obtained from the inhomogeneous density profile of the trapped atomic cloud. The 
wings of the cloud realize the high field limit where the polarization approaches f00%, which is not 
accessible for an electron gas. 
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I. INTRODUCTION 

Ultracold atoms can be prepared with almost complete 
control over their density, temperature and interactions. 
They serve as model systems for exploring unsolved prob- 
lems in many-body physics as well as for demonstrat- 
ing well-known textbook physics, such as of ideal, non- 
interacting gases that do not exist in nature. For example 
several phenomena related to Pauli blocking of fermions 
were clearly demonstrated only after the advent of ul- 
tracold Fermi gases Here we use the tunability of 
atomic interactions near Feshbach resonances to create a 
non-interacting Fermi gas with two components. This is 
a realization of an ideal, non-interacting free electron gas 
(FEG) with spin up and down components, as assumed 
in the simple theory of metals. We demonstrate how 
such an ideal Fermi gas will respond to effective mag- 
netic fields, which is described by Pauli paramagnetism. 

The paper is mainly pedagogical. It explains how para- 
magnetism is observed in trapped atomic samples which 
have an inhomogeneous density due to the harmonic con- 
finement potential. Furthermore, these atomic samples 
realize a canonical ensemble in the fixed atom numbers 
iVf and of the two components, whereas in metals 
and solid-state physics, Nf — JVj, is the magnetization 
determined by the applied external magnetic field. In 
all previous studies of paramagnetism, the magnetiza- 
tion was weak since the applied magnetic field times the 
magnetic moment was much smaller than the Fermi en- 
ergy With ultracold atoms, we can easily realize 
the strong field case where the chemical potential differ- 
ence is larger than the Fermi energy and therefore fully 
polarizes the gas. Besides its pedagogical purpose, this 
paper experimentally demonstrates Pauli paramagnetism 
in a truly non-interacting system exactly described by ba- 
sic theory, whereas measurements even in simple metals 
revealed major discrepancies due to interaction effects 
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FIG. 1. Comparison between a FEG and a trapped two- 
component Fermi gas at temperature T = 0. (a) Response 
of the FEG to an external magnetic field B. Electrons in 
the two spin states are filled up to the Fermi energy Ef- An 
external magnetic field B shifts the energy of the two states 
into opposite directions and produces the density difference, 
or magnetization, (b) Density distributions n-f and nj_ of a 
two-component trapped Fermi gas at T = 0. As discussed in 
the text, the local density polarization and corresponding nor- 
malized effective magnetic field field varies across the atomic 
cloud. This is illustrated in the diagrams next to the curve 
for n-j-, where the gaps between the two spin states represent 
the absence of spin-flip processes in the system mentioned in 
the text. 



II. CONCEPT 
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Pauli paramagnetism explains the magnetization of a 
free electron gas with two spin states, neglecting the con- 
tribution from orbital motion. In a FEG, the applied field 
shifts the energies of the two spin states in opposite direc- 
tions as in Fig.fJJa). For a uniform FEG at temperature 
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FIG. 2. Determination of the susceptibility \ from atomic 
density profiles, (a) Simulated density profiles of a two- 
component Fermi gas in a spherically symmetric harmonic 
trap at zero temperature. The total population imbalance 
5 = (JV t - N.i)/{N t + Ni) = 21%. The shaded area cor- 
responds to the partially polarized region where x is well- 
defined, (b) The local chemical potential difference Afi with 
(solid) and without (dashed) normalization by the Fermi en- 
ergy Ep[n]. For the fully polarized gas A/i = A/j,/Ef[ti] = 1 
according to the definition in the text, (c) The density dif- 
ference An with (solid) and without (dashed) normalization 
by the total density n. (d) Normalized density difference An 
vs. normalized chemical potential difference Afi. For the 
fully polarized gas, An = A/i = 1. A linear fit at low A/i 
(dashed line) gives \ = 3/2 4/3 « 1.19. (e) Normalized mag- 
netic susceptibility x versus the normalized chemical potential 
difference Afi. 



T=0, the total number N^m of electrons in each spin 
state inside a volume V is given by 



N tU) = 



g(E ± u B B)dE, 



(1) 



where g{E) = fa ) 3/2 E 1 ' 2 is the density of states, 
Ep the Fermi energy, m the mass of an electron, n B the 
Bohr magneton and h Planck's constant. The magnetic 
susceptibility x 1S defined as usdAn/dB, where Art = 
rif — nj. = (N^ — N±)/V. As a dimensionless quantity, we 
consider the normalized susceptibility 



X 



d (An/n) 



fi B d(B/E F [n/2] 



(2) 



where the polarization An/n is used as a more meaning- 
ful measure of magnetization. Ep[n] = h 2 (6TT 2 n) 2 ^ 3 /2m 
and the total density is n = + n± = (Nf + N±)/V. 

Here we experimetanlly simulate Pauli paramagnetism 
using trapped fermionic alkali atoms. The two low- 
est hyperfme states of 6 Li, \F = 1/2, m F = 1/2) = |f) 
and |F = 1/2, mp = —1/2) = ||), represent the two spin 
states of an electron. This atomic system is different 
from free electron systems: (1) due to the slow dipolar 
relaxation the population in each spin state is con- 
served and therefore, (2) an external magnetic field does 
not lead to any magnetization, it only shifts the energies 
of the two spin state (only in the case of strongly mag- 
netic dipolar gases such as Cr the dipolar relaxation is 
fast, and therefore the spin populations follow an exter- 
nal magnetic field @). 

To realize the effective magnetic field we introduce den- 
sity imbalance while preparing the system 0, H| as illus- 
trated in Fig. Hlb). I n the grand-canonical description 
of a system with fixed atom number, a Lagrange mul- 
tiplier (called the chemical potential) ensures the cor- 
rect expectation value for the atom number. In our case, 
the atom numbers in both spin states are fixed. In a 
grand-canonical description, this leads to two Lagrange 
multipliers (the two chemical potentials for the two spin 
states), or to one Lagrange multiplier for the total num- 
ber, and one for the population difference. It is the sum of 
the latter Lagrange multiplier and the external magnetic 
field which constitutes the effective magnetic field in our 
ensemble. The contribution from the external magnetic 
field is always canceled out by an additive component 
in this multiplier. For example, a balanced gas always 
has zero effective magnetic field independent of the ap- 
plied field. Therefore the effective magnetic field depends 
only on the population imbalance and not on the external 
magnetic field. 

For a uniform Fermi gas at T = 0, the number in each 
spin state determines chemical potential, A«-f(n @ and 
satisfies 



t(-l-) 



g(E)dE= / g(E ± Au/2)dE, (3) 



TUbB 



where fi = (fif + ni)/2 and Afi = fif — a^. g(E) is the 
same as in Eq. (TTJ) if m is taken to be the mass of a 6 Li 
atom. Comparing with Eq.(fT]) shows that the chemical 
potential difference A/i is the effective magnetic field that 
takes the place of the magnetic field in a FEG. With this 
effective magnetic field, the magnetic susceptibility can 
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be written as 



III. EXPERIMENT AND RESULTS 
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We show now how a single atomic density profile in a 
harmonic trap V(r) can be used to determine the nor- 
malized susceptibility as a function of magnetic field. 
In the local density approximation, the density at each 
point in the trap is that of the corresponding uniform gas 
with modified chemical potentials /i^)(r) = M^m — V(r) 
0) EH > where m are global chemical potentials defined 
for the whole trapped cloud which constrains the total 
number of atoms in each spin state. This concept has 
been used to determine the equation of state of interact- 
ing Fermi gases [H, [111-1131 . Here we are interested in the 
normalized magnetic susceptibility which is the slope of 
the normalized density difference versus the normalized 
chemical potential difference [3]. We demonstrate this 
procedure first using simulated density profiles (Fig. [5]). 

As ftfd) (f) varies across the cloud, the chemical po- 
tential difference A/i(r) = /if(r) — = /i£ — /if re- 
mains constant, as shown in Fig. [UJb). However, the 
total density n(r) = n^(r) + n±(r) changes across the 
cloud, and consequently the local Fermi energy Ep[n(r)] 
as well, with -EfM being the same as in Eq. (TTJ) if m is 
taken to be the mass of 6 Li. The normalized magnetic 
susceptibility % depends on the dimensionless effective 
magnetic field A/i = A/i(f) / Ep[n(r)], which does vary 
across the cloud, from small values near the center to the 
saturated value of one at the edges. From Fig. [21c), a 
single set of atomic density profiles contains a range of 
polarizations An = An/n and normalized effective mag- 
netic fields Afj,/Ep[n]. The slope of the An vs. An yields 
the normalized susceptibility 



X = d(An) jd (Am) 



(5) 



The normalization of \ in Eq. ([5]) is chosen to ensure 
An = A/i = 1 for a fully polarized gas. The slope near 
the origin, corresponding to small values of A/i, gives 
X(T = 0) = x = 3/2 4 / 3 w 1.19 [15]. As shown in 
Fig. [He), the susceptibility decreases with increasing A/i 
and drops to zero when the Fermi gas reaches full polar- 
ization. 
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FIG. 3. Double-shot phase contrast imaging. |e) is the 2P3/2 
excited state in 6 Li. Two successive images are taken at two 
different imaging frequencies, (a) and (b). Corresponding 
CCD images are shown on the right. 



An ultracold one-component Fermi gas of 6 Li in the 
state \F = 3/2, mp = 3/2) is prepared by sympathetic 
cooling with bosonic 23 Na atoms [l6|,[I3|- 6 Li atoms are 
then loaded into an optical dipole trap followed by radio- 
frequency (rf) transfer to \F = 1/2, mp — 1/2} using 
an rf Landau-Zener sweep at 300 G. A superposition of 
the two lowest hyperfine ground states, \F = l/2,mp = 
1/2) = |t) and \F = 1/2, m F = -1/2) = \\) is prepared 
by a rf sweep at 300G, and holding here for 500 ms results 
in an incoherent mixture. Further cooling is provided by 
lowering the optical trapping potential which leads to 
radial evaporation (axial confinement is provided mag- 
netically). We increase the magnetic field in 200 ms to 
528G, where the scattering length is zero [l8j], leading to 
a non-interacting two-component Fermi gas. The oscilla- 
tion frequencies of the final trapping potential are u) y — 
2tt x 35 Hz axially and lo x = lj z = 2ir x 390 Hz radially. 
The total population imbalance S = (Nf-N i )/(Nf+N i ,) 
is about 33%. The final temperature is about 0.22T/Tp 
where ksTp — Ep[n/2] with n being the total peak den- 
sity and ks the Boltzmann constant. 

To obtain density profiles of each spin state of 6 Li, 
we use a double-shot phase contrast imaging technique 
that involves two images taken at two different imaging 
frequencies d, Oj} (see Fig. [3]) . One frequency is tuned 
between the two resonant transition frequencies for \\) 
and \l) states to the 2P 3 / 2 electronically excited state. 
The resulting phase contrast picture is a measure of the 
density difference An = n-j- — nj,. To measure the density 
of each component, a second image is taken of the same 
cloud at a frequency detuned towards the \f) transition 
by 10 MHz. The two images can be taken in quick suc- 
cession because of the non-destructive nature of phase 
contrast imaging. 

For dispersive imaging techniques such as phase con- 
trast imaging, it is crucial to carefully focus the imaging 
system, since the atomic cloud refracts the probe light (in 
contrast to resonant absorption imaging). One can check 
this by comparing in-trap profiles of |t) atoms imaged at 
positive and negative detunings of 60 MHz and find the 
focal position which minimizes the difference [l7| . 

Figure 2] shows the experimental results. The camera 
images provide two-dimensional (2D) column densities 
712D which integrate the density distribution along the 
line of sight. To reduce the noise for the following analysis 
we take advantage of the symmetry of the trap and per- 
form quadrant averaging and elliptical averaging: First, 
the four quadrants around the center of the images are 
averaged. Second, the data is averaged along elliptical 
contours. The aspect ratio of the ellipses is determined 
from in-trap images [§[. Temperature and global chemi- 
cal potential for each component are obtained by fitting 
2D densities with theoretical non-interacting density pro- 
files. Since the temperature T is small but non-negligible, 
data are fitted with finite-temperature Fermi gas 2D dis- 
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FIG. 4. (color online) Experimental determination of x- Rep- 
resentative profiles of (a) 2D column densities for |f) (upper 
solid curve) and |4_) (lower solid curve). Temperature and 
chemical potential are obtained from the fit (dotted curves) 
as discussed in the text. Each elliptical average is plotted ver- 
sus the axial position (long axis of the ellipse) . (b) 3D density 
difference obtained from a single double-shot image. Dotted 
curve is a theoretical curve based on the chemical potentials 
and the temperature obtained from the 2D fit in (a), (c) The 
normalized 3D density difference is plotted as a function of 
the normalized chemical potential difference. The solid line is 
the theoretical prediction, (d) The linear fit at low A/i gives 
the low-field susceptibility x=0.95(l). 



along the line of sight. Since the majority component has 
better signal-to-noise ratio, we determine the tempera- 
ture and chemical potential first for this component, and 
keep the temperature as a fixed parameter in the fit of 
the minority profiles. Three-dimensional (3D) densities 
[Fig. SJb)] are obtained by performing the inverse Abel 
transformation of the column densities Q . The linear fit 
of An as a function of A/z at small A/z gives the suscep- 
tibility at finite temperature [Fig. 0|d)]. For finite but 
low temperatures, the susceptibility is known to vary as 



m 



X = Xo 



1 - 



12 \ T, 



T 



(7) 



T/Tp changes from 0.24 to 0.35 within the fitting re- 
gion in Fig. Hfd), affecting the susceptibility by about 
3%. For T/Tp=0.29, this theoretical calculation gives 
X = 1.108 and the experimentally obtained value is 
0.95(1) in Fig.@|d). The discrepancy between the experi- 
mental and theoretical value is most likely due to a resid- 
ual dispersion effect lea ding to systematic uncertainty es- 
timated to be 10-20% Ha [HI- We did not attempt to 
quantify this effect since this paper is mainly pedagogi- 
cal and a dispersion effect depending on the imbalance as 
well as two detunings would be time-consuming to quan- 
tify 

Figure U{c) illustrates that we can easily approach the 
high field limit of a fully polarized gas A[i w 1. Note that 
An is smaller than 1 because at finite temperature the 
minority component extends beyond its Thomas-Fermi 
radius. In a metallic free electron gas, with a typical 
Fermi temperatures of 10,000 K, the high field region 
would require magnetic fields of 10,000 Tcsla, 200 times 
stronger than the strongest continuous laboratory mag- 
netic fields. 



IV. CONCLUSIONS 

We have demonstrated how the ideal gas magnetic sus- 
ceptibility can be measured in ultracold Fermi gases. Af- 
ter considering finite-temperature effects, the experimen- 
tal results are in quantitative agreement with prediction 
from Pauli paramagnetism. Since the local spin polariza- 
tion varies across the trapped atomic cloud, a single den- 
sity profile gives the susceptibility for both low and high 
effective magnetic fields. This method can be applied to 
the study of magnetic properties of strongly interacting 
Fermi gases which are of current interest [22l425| . 



tribution in a harmonic trap given by 

n 2 u(x,y) = 
m(k B Tf 



-Li 2 



-rn(cj 2 n:X 2 +u 1 2 y v 2 )/2)/kBT) 



> (6) 



27th 3 LU z 

where Li„(z) is the nth-order Polylogarithm. z is the axis 
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